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Abstract 



Directed sideward flow of light charged particles and intermediate mass frag- 
ments was measured in different symmetric reactions at bombarding energies 
from 90 to 800 AMeV. The flow parameter is found to increase with the charge 
of the detected fragment up to Z = 3 — 4 and then turns into saturation for 
heavier fragments. Guided by simple simulations of an anisotropic expand- 
ing thermal source, we show that the value at saturation can provide a good 
estimate of the flow angle, ©fl ow , in the participant region. It is found that 
0fiow depends strongly on the impact parameter. The excitation function of 
©flow reveals striking deviations from the ideal hydrodynamical scaling. The 
data exhibit a steep rise of ©fl ow to a maximum at around 250 — 400 AMeV, 
followed by a moderate decrease as the bombarding energy increases further. 

Keywords : Heavy-ion collisions, Reaction plane, Directed sideward flow, Flow parameter, 
Flow angle, Expanding thermal source, hydrodynamical scaling. 
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I. INTRODUCTION 



The emergence and development of a collective expansion motion in the course of an 
energetic heavy-ion collision reflects the response of the nuclear system to the internal pres- 
sure built-up at the interface of the interacting nuclei. Thus, experimental investigations of 
this phenomenon are expected to provide information about the compressibility coefficient 
characterizing the nuclear equation-of-state which is one of the major objectives in nuclear 
physics research. This has generated a widespread interest in this subject which is currently 
the object of intensive experimental activities [|TJ over a broad range of bombarding ener- 
gies going from the region where a possible liquid-gas phase transition might occur to the 
ultra-relativistic energy domain where the transition to the quark gluon plasma is expected 
to take place. 

In finite impact parameter collisions, because of the presence of cold spectator remnants 
the expanding nuclear matter in the hot and dense central region is deflected sidewards 
along a preferential emission direction in the reaction plane. The most natural experimental 
observable to characterize this effect is the flow angle, Ofl ow , &-e. the angle between the 
direction of the collective sideward motion and the longitudinal beam axis. This observable is 
of relevance not only because of its sensitivity to the stiffness of the nuclear matter equation- 
of-state, but its knowledge is also essential in other measurements such as the extraction 
of the squeeze-out signal around the flow axis |||-§J and two-particle interferometry ||. Its 
determination requires the use of "47r" detectors capable of global event reconstruction. 

The sideward flow can be characterized using the transverse momentum analysis 
method ||. The latter method allows one to extract the so-called "flow parameter", F$, 
which is defined as the slope at mid-rapidity of the average in-plane transverse momentum 
as a function of the rapidity. This observable is by definition representative of the highly 
excited and compressed central region. It has also the advantage that uncertainties on the 
reaction plane determination due to finite number of particle effects and detector biases 
can be accounted for by introducing an appropriate correction factor |6],|7|]. An inherent 
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problem in this method is that the flow parameter does not reflect only the collective be- 
haviour but it is also sensitive to the random thermal motion of the emitted particles. In 
the absence of thermal fluctuations the scaled flow parameter = Fg/p^ m -, where p^ m ' is 
the projectile center-of-mass (cm.) momentum per nucleon, is a measure of the flow angle 
(F s (0) = tan(e flow )), while in the presence of a large thermal motion the value of might 
be significantly lower than tan(Bfj ow ). It follows thereby that intermediate mass fragments 
(IMFs), which are less subjected to thermal fluctuations than do lighter particles, are ex- 
pected to be strongly aligned along the flow direction. This effect was first observed by 
the Plastic Ball group || and confirmed later on in several experiments P,|9|-[13|. Recently, 
quantitative measurements of IMF sideward flow in the Kr+Au reaction at 200 AMeV 
have shown that the flow parameter increases with the size of the detected particle and 
reaches a constant limiting value for fragments with masses 4 < A < 12 These con- 
siderations underscore the importance of accurate measurements of the flow parameter of 
IMFs, in particular at beam energies of a few hundred AMeV where IMFs are copiously 
produced |D|JT|-fB|. 

Many interesting processes might simultaneously contribute to the observed flow : the 
release of compressional energy [I7|, the thermal pressure ITS], the momentum dependence of 
the nuclear force [Jl9"| , p0f and the in-medium nucleon-nucleon cross sections |2]J . Systematic 
studies by varying the bombarding energy and the system size are needed to disentangle 
the different contributions of these effects since they are expected to exhibit different depen- 
dences on the initial conditions. 

In the present paper we report on detailed experimental results of the flow parameter of 
light charged particles and IMFs measured in different symmetric reactions, Ru+Ru, Xe+Csl 
and Au+Au, at beam energies between 90 AMeV and 800 AMeV. We show that an accurate 
estimate of the mean flow angle can be obtained from the observed flow parameter of heavy 
IMFs. Then, by applying this idea we investigate the dependences of the flow angle as a 
function of the collision impact parameter and the bombarding energy. 
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II. EXPERIMENTAL SETUP 



The experimental results reported in this paper were obtained in a series of experiments 
carried out using the FOPI detector [^2J at the SIS accelerator facility of GSI-Darmstadt. 
Several symmetric reactions at beam energies going from 90 AMeV to 1.9 AGeV were 
investigated. Here, we will present data for the following systems : Ru on Ru and Xe on Csl 
at E = 400 AMeV and Au on Au at E = 90, 100, 120, 150, 250, 400, 600 and 800 AMeV. The 
beam intensities were typically 10 5 ions/s. The target thickness was between 100 mg/cm 2 
(at the lowest incident energy) and 400 mg/cm 2 (at the highest energy), corresponding to 
an interaction length going from 0.5% to 2%. Data at 600 and 800 AMeV were obtained 
from the first generation of FOPI experiments using the Phase I setup, while the lower 
beam energy data were taken more recently with the Phase II configuration of the FOPI 
detector. Details about the FOPI apparatus and its performances have been reported in 
previous publications [0,[23[]. Here, we recall briefly some of the sub-detector components 
of particular interest in the present analysis. 



In its Phase I configuration |[22|| , the FOPI detector consisted of a highly segmented 
Forward Wall of 764 plastic scintillators, divided into an external (512 strips) and an internal 
(252 scintillator paddles of trapezoidal shape) components. To identify the slow heavy 
fragments stopped in the external wall an ensemble of 16 large gas ionisation chambers with 
a total of 128 anodes of individual readout, positioned in front of the external wall, was also 
used. For the same purpose the inner part was supplemented with a shell of 60 thin (2 mm) 
plastic scintillator paddles. The whole setup covered the laboratory polar angles, ©i a b, from 
1.2° to 30° over the full azimuth. This device allowed us, event by event, to identify the 
nuclear charge and measure the vector velocities of most of the light charged particles and 



IMFs (up to Z = 12) emitted in the forward cm. hemisphere. In phase II [^3| experiments, 
a cylindrical drift chamber CDC mounted inside a superconducting solenoid was used at 
backward angles (©i a b = 30° - 150°). Pions, protons and deuterons were identified in the 
CDC by means of their mean energy loss < dE/dx > and their laboratory momentum, 
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obtained from the curvature of the particle tracks in the field of a 0.6 T strength. The high 
granularity of the setup allowed high multiplicity events to be measured with a negligible 
multi-hit rate. The apparatus ensured also a very good azimuthal symmetry which is an 
important feature for the study of the flow phenomenon. As the main objective in this work 
was to extract information on the sideward flow of IMFs, our analysis was based mainly on 
the data from the forward sub-detectors. The CDC was used only for the purpose of event 
characterization, as we will see in the next section. 

All data presented in this paper were obtained from the analysis of events taken under 
the "central trigger" condition ||22|| . The latter was defined by adjusting the charged particle 
multiplicity to a value which corresponds to impact parameters less than ~ 2/3 of the 
maximum impact parameter. At each bombarding energy, samples of about 10 5 to 10 6 of 
such events were recorded. This large amount of available events allowed us to extract high 
statistics data in particular for IMFs. 

III. EVENT CHARACTERISATION 

To extract quantitative information on flow phenomena from the data, one requires a 
good event characterization both in impact parameter and azimuth of the reaction plane. 

In order to classify the measured events according to their degree of centrality, we have 
employed the standard method based on the correlation between the multiplicity of the 
emitted particles and the impact parameter. The event multiplicity was extracted as the 
number of charged particles detected per event in both the outer part (from 7° to 30°) of 
the forward scintillator wall and the CDC. For Au+Au at 600 and 800 AMeV where the 
CDC was not operational, the event multiplicity was restricted only to charged particles 
measured by the outer part of the scintillator wall. In the following, the event multiplicity 
will be labelled MUL when the CDC is included and PMUL otherwise. The measured event 
multiplicity distributions exhibit the typical plateau for intermediate values followed by a 
rapid fall off for the highest multiplicities |T^,|T^[. These distributions were divided into five 
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intervals, in accordance to the procedure introduced in previous works fl5| , |24]| . The highest 
multiplicity bin, named MUL5 (or PMUL5) starts at half the plateau value and the remaining 
multiplicity range was divided into four equally spaced intervals, named MUL1 to MUL4 (or 
PMUL1 to PMUL4). Most of the results which will be presented in this paper deal with the 
MUL4 (or PMUL4) event class corresponding to an average geometrical impact parameter of 
~ 3.5 fm (for the Au+Au reaction). As we will see, this centrality bin lies in the region where 
the sideward flow reaches its maximum (|T~3| , [T4|1 . The corresponding cross sections and mean 
geometrical impact parameters (obtained by assuming a sharp-cut-off approximation) are 
given in Tab. | and Tab. [n]. Note in passing that the value of the reduced impact parameter 
is nearly the same in all cases which will allow direct comparisons of the experimental results 
obtained for different systems and at different bombarding energies. 

To reconstruct the reaction plane, we have used the transverse momentum analysis 
method ||. In order to remove autocorrelation effects, the azimuth of the reaction plane was 
estimated for each particle i in a given event as the plane containing the vector Q 1 and the 
beam axis where Q 1 is calculated from the transverse momenta p[ of all detected particles 
except the particle of interest i 



M 



Q' 1 = $>Vt+mH)- (1) 



j=i 



M is the multiplicity of the event and cj j = 1 if y^ > 5, — 1 if y^ < —5 and otherwise, 
is the j th particle rapidity divided by the projectile rapidity in the cm. system. The 
parameter 5, choosen equal to 0.5, was introduced in order to remove the contributions of 
mid-rapidity particles which have a negligible correlation with the reaction plane. According 
to Ref. ||25|| , a boost velocity v l h = p\/(rrf ys — m 1 ) (m 1 is the mass of particle i and m sys is 
the sum of the projectile and target masses) was applied to each particle j in order to take 
into account the effects of momentum conservation due to the exclusion of the particle of 
interest i. In order to estimate the accuracy on the reaction plane determination, due to 
finite number of particle effects and detector biases, we have used the method described in 
reference || which consists in randomly subdividing each event into two and calculating on 



average the half difference in azimuth, A$r, between the reaction planes extracted from 
the two sub-events. A$r gives an estimate of the dispersion of the reconstructed reaction 
plane with respect to the true one ||. The results are displayed in Tab. | and Tab. [TI| in 
terms of the standard deviation width ct(A$r) extracted from a gaussian fit to the A$r 
distributions. As can be seen, the reaction plane is, in all cases, rather well estimated, with 
a precision which varies typically from ~ 22° to ~ 46° depending upon the system and 
the incident energy. All data subsequently presented in this paper are corrected for these 
uncertainties on the reaction plane determination. 

IV. EXPERIMENTAL RESULTS AND DISCUSSION 

In the framework of the transverse momentum method || , the sideward flow is quantified 
by plotting the mean in-plane transverse momentum per nucleon < p x > as a function of 
the rapidity. For each particle, the quantity < p x > is obtained by projecting its transverse 
momentum p[ onto the reaction plane (the x direction is defined to be in the reaction 
plane). The reaction plane is reconstructed by removing autocorrelation effects as described 
in the previous section. In the present work, the data have been expressed in terms of 
scale invariant (dimensionless) quantities in order to avoid trivial scaling with the incident 
beam energy [p6f| . In what follows, scaled quantities will be indicated by the index (°/ : 
< p x ' > = < ^ > will denote the mean in-plane transverse momentum per nucleon < p x > 

Pp 

scaled to the projectile momentum per nucleon p 1 ^ 1 in the cm. system and = will 
refer to the cm. particle rapidity y cm normalized to the cm. rapidity yT 1 of the system. 

Displayed in Fig. 1 are representative examples of < p x ^ > versus y^ plots for the three 
reactions Au+Au, Xe+Csl and Ru+Ru, all at a beam energy of E = 400 AMeV. The 
data belong to the multiplicity bin MUL4 where, as outlined earlier, the directed sideward 
flow is close to its maximum. The < p x ^ > values have been divided by the mean cosine of 
A$r in order to take into account the resolution of the reaction plane determination. This 
has been done according to the method recently proposed by Ollitrault [0. The values of 
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< cos(A$r) > are given in Tab. | and Tab. |TJ. The small statistical uncertainties (smaller 
than the symbol sizes for Z = 1 and Z = 2 particles) on the experimental points, in particular 
for IMFs, illustrate the high statistics recorded in the present work. For the Ru+Ru reaction, 
the data for fragments with Z > 4 are subjected to relatively large statistical errors. There 
are two reasons for that : i) for this system only half of the available statistics was used and ii) 
the production yields of IMFs are lower in lighter systems. The plots of Fig. 1 were obtained 
from the analysis of the data taken by the sub-detector components composing the forward 
Wall (see section II), allowing IMF measurements. This is the reason why the plots of Fig. 1 
are shown only in the forward hemisphere in the cm. frame (positive cm. rapidities). It is 
worth recalling that the forward Wall provides an individual element identification of light 
charged particles and IMFs (up to Z = 12) along with velocity measurements. Thus, the 

< p^ > versus distribution measured for a given fragment includes the contributions of 
all associated isotopes. It should be, however, stressed that the lack of mass identification 
does not affect the < p^ > quantity as it is directly derived from the measured velocities 
without any assumption on the mass of the particle. 

As shown in Fig. 1, the measured < p^ > versus distributions exhibit the typical 
S-shape behavior ||^7| reflecting the transfer of momentum between the backward and 
forward hemispheres. The linear part of the S-shaped curve is representative of the flow 
in the participant region, the so-called side-splash effect, while the fall-off observed at high 
rapidities is caused by the bounce-off effect |p8| , p9|l . It is customary to quantify the magnitude 
of the participant flow as the slope of the < > versus y^ curve at mid- rapidity [^7]] 



Ff = d<pf >/dy (o)|^ (2) 

In practice, we extracted this quantity (known as "flow parameter" in the literature) 
by fitting a polynomial function of the form a + x y^ + c x (y^) 3 to the data. The fit 
was restricted to the linear branch of the S-shaped curve. We have verified that reasonable 
changes of the limits of the y^ range where the fit is applied, lead to the same values of 
.Fg '* within the statistical error bars. It is worth to notice that the influence of momentum 



conservation effects on the flow parameter was found to be very weak (less than a few %). 
More details on the fitting procedure can be found in Ref. ||f§]. 

Before presenting the results, let us briefly comment on the influence of the biases intro- 
duced by our apparatus on the measured values of the flow parameter. To evaluate these 
effects, we performed detailed simulations where theoretical events were passed through the 
FOPI detection filter including geometrical cuts and energy thresholds. The theoretical 
events were obtained from Quantum Molecular Dynamics (QMD) model calculations using 
the so-called IQMD version H30| , |3i1| which yields generally a quite good agreement with the 
sideward flow data |3^j33|] . We estimated the effects of detector cuts by comparing the flow 



parameter calculated with and without including the experimental filter. Because of limited 
statistics (we used samples of a few hundred IQMD events per impact parameter unit), 
this could be done accurately only for light particles. We found that the observable is 
mainly affected by the ©i a b = 30° cut. The results indicate that acceptance effects depend 
on the fragment charge, the collision centrality and the beam energy. For MUL4 events at 
an incident energy of E = 250 AMeV, was found to be biased down by about 30% and 
10% for Z = 1 and Z = 2 particles, respectively. For heavier fragments, these effects are 
expected to be much lower as the 0i a b = 30° cut has a weaker influence on IMFs. This 
will be illustrated later on (Fig. 5) using more simple simulations based on the decay of an 
expanding thermal source. 

A. Centrality dependence of the flow parameter 

Figure 2 shows the dependence of the scaled flow parameter measured for 
Au(400 AMeV)+Au as a function of the collision impact parameter. Error bars corre- 
spond to statistical uncertainties multiplied by to take into account the uncertainty 
from the polynomial fit to data. The impact parameter was obtained from the measured 



cross sections (b geo = ycr/V) assuming a sharp-cut-off approximation. The flow parameter 
F^ is presented here in the form of a coalescence invariant quantity |3"B|-|3~5] , i.e. including 
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the contributions of all detected particles each weighted by its measured charge. 

The results exhibit a maximum of the sideward flow around b = 4.8 fm. The observed 
trend is qualitatively consistent with the earlier results [jSTJ reported by the Plastic Ball 
group at the same incident energy. It can be intuitively understood from the expectation 
that the flow parameter should be zero at b = 0, for symmetry reasons, and must tend also 
towards zero in peripheral collisions. The value of at maximum (~ 0.38) is, within 
acceptance effects, in good agreement with the data of the Plastic Ball group and those 



more recently published by the EOS Collaboration ]10 



B. Fragment charge dependence of the flow parameter 

The fragment charge dependence of the flow parameter is depicted in Fig. 3 for Au+Au 
reactions at different bombarding energies. Error bars correspond to statistical uncertain- 
ties multiplied by \px 1 . The data are not corrected for the distorsions introduced by the 
FOPI apparatus. The influence of detector cuts on the measured flow parameters has been 
already discussed. Qualitatively, the data exhibit the same pattern at all five beam energies, 
characterized by a gradual increase of with the charge of the detected particle followed 
by a clear tendency to level-off above a certain value of Z. With increasing beam energies, 
the satuation appears at lower Z values. This is illustrated by the curves representing the 
results of a fit with a Fermi function to the data. 

A similar behaviour has been recently reported in the Kr + Au reaction at 200 AMeV 
where the flow parameter was found to reach a constant limiting value for fragments with 



masses 4 < A < 12 [11]. It is worth-while mentioning that this trend has been also seen for 



the squeeze-out effect P,|36[| and seems to be typical for all flow observables ]13 



The increasing strength of the sideward flow with the fragment size was observed for 
the first time || by the Plastic Ball group for Au(200 AMeV)+Au reactions and confirmed, 
since then, in several experiments p|,p|-4T^1. This phenomenon was predicted by hydrody- 



namical models 13711 where collective effects were found to be much more visible for heavier 
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mass fragments. QMD calculations also yield larger flow for heavier fragments ||38|| . This 
large flow carried by fragments is of particular interest because of its enhanced sensitivity 
to the parametrisation of the nuclear equation-of-state used in dynamical model calcula- 



tions PLpBp) 



It can be also understood from the interplay between collective and thermal (random) 
motions. In an idealized picture in which nucleons and fragments are emitted from a common 
expanding thermalized source, heavy fragments are less subjected to thermal fluctuations 
than do lighter particles. Thermal fluctuations are governed by the thermal energy which 
is independent of the mass of the particle, while the collective expansion energy increases 
linearly with the mass and is, therefore, better reflected in heavier fragments. Thus, it is 
easy to understand, within this simple picture, that in the presence of thermal fluctuations 
the apparent flow angle that one can extract from the flow parameter 0^ = arctan(Fs < -°' ) ) is 
lower than the effective flow angle. The nearly constant value observed in the flow parameter 
of heavy fragments (Fig. 3) can be attributed to the fact that above a certain mass (or charge) 
the emitted fragments are only very little affected by the random thermal motion so that 
their apparent flow angles are expected to be very close to the effective flow angle. That is 
the idea that we will exploit below (section IV. D) in order to extract an estimate of the flow 
angle from the data. 

It should be mentioned that it was shown in Ref. |34j] that the increase of the sideward 
flow with fragment size can be also described by a simple momentum space coalescence 
prescription providing a transverse momentum cut of 0.2 AGeV is imposed. 



C. System size dependence of the flow parameter 

In Fig. 4, the normalised flow parameter as a function of the charge of the detected 
fragment is displayed for the three reactions under study. The flow parameter is divided here 

1/3 1/3 

by the quantity A P + A T , where Ap and At denote the mass of, respectively, the projectile 
and target nuclei. Two observations can be readily made from the examination of this figure. 
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First, the dependence of the flow parameter on the fragment charge measured in Ru+Ru and 
Xe+Csl reactions follows the same trend as observed in the heavier system Au+Au. Second, 
the projectile-target mass dependence of is consistent, within statistical uncertainties, 
with the Ap 3 + A^ 3 empirical scaling rule introduced recently by the EOS collaboration [^0 . 
The origin of this scaling rule is not yet well understood. It might be attributed, in a 
hydrodynamical picture, to the fact that for collisions with velocities well above that of 



sound, the pressure built-up should scale with collision length or time [(40 



D. Extraction of the flow angle from IMF measurements 

As outlined in the introduction, the most natural observable characterizing the sideward 
deflection of the nuclear matter emitted in non-zero impact parameter collisions is the flow 
angle i.e. the angle between the direction of the collective sideward motion and the longitu- 
dinal beam axis. Indeed in contrast to the flow parameter, 0fl ow is a measure of the overall 
emission direction of all particles belonging to a given event and is not affected by thermal 
fluctuations. It follows therefore that the determination of this observable, because of its 
greater sensitivity, should help to better understand the origin of the flow and disentangle 
the different phenomena that may contribute to the observed effect. The knowledge of 0fl ow 
is, on the other hand, of great importance in the investigation of other phenomena such as 
the out-of-plane squeeze-out effect [@-|4j and two-particle interferometry . 

Several methods were introduced to reconstruct this observable in high energy heavy-ion 



experiments. The sphericity method [|41j| , based on the diagonalization of the momentum 
flow tensor, allows the extraction of the flow angle as well as two aspects ratios characterizing 
the event shape (assumed to be ellipsoidal). This procedure has the advantage that it pro- 
vides an event-by-event shape characterization. It deals, however, with the overall emission 
pattern including the spectator component. Furthermore, the sphericity analysis is strongly 
affected by the distorsions due to the effects of finite number of particles ||2] . An alternative 
shape analysis method, where the contribution of the spectator matter can be removed, was 
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proposed by Gosset et al |43[]. The flow angle is adjusted, within this method, by fitting a 



simple anisotropic gaussian distribution to the triple differential momentum distributions of 
particles detected in the participant region. This procedure is nevertheless reliable only in 
the case of low impact parameter collisions where the contribution of the spectator matter 
is not very important jPJ . 

We propose in the present article a new method which exploits the saturation observed 
in the dependence of the flow parameter as a function of the fragment charge (Fig. 3 and 
4). As discussed earlier this saturation is due to the fact that IMFs are much more aligned 
along the flow direction than do lighter particles : "IMFs go with the flow" |39],fi5|. One 
expects therefore that if the phase space region occupied by IMFs is sufficiently elongated 
then their apparent flow angle 0^ should give a good measure of the effective flow angle. 

To investigate more quantitatively this idea, we have performed Monte-Carlo simulations 
based on an anisotropically expanding thermal source calculations, where particles with 
different masses share the same thermal energy but their collective energy is proportional 
to their mass. We assume that, at the freeze-out stage of the collision, the nuclear system 
is in local thermal equilibrium, i.e. same temperature throughout the entire volume of the 
source. The latter is considered as a cylinder, in configuration space, whose principal axis 
coincides with the flow axis. The flow angle is introduced as a free parameter, Og^ 1 *, in 
the simulations. Particles are generated one by one independently of the influence on each 
other. The momentum of a given fragment is considered as resulting from the superposition 
of a purely collective component (common velocity field to all fragments) on the top of a 
thermal motionP] 

p = p th + p co11 (3) 
The random thermal component p* h is assumed to obey a Maxwellian distribution 

M(p th ,m) = p th V^ m2+ P th2 / T (4) 



L details about the simulations can be found in Ref. 
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where m is the mass of the fragment and T is the temperature of the source. 

The collective component p coYl is calculated assuming that at the moment of the explo- 
sion, the expansion velocity of a fragment increases linearly along the spatial coordinates 
towards the surface of the freeze-out volume 



= n 2 < EgVA > 
q, V 931.8 I 

J i=x,y,z 

The fragment positions r ; (i = x, y, z), expressed in the source reference frame, are 
randomly determined inside the cylindrical volume of the source assuming a uniform density 
at the freeze-out. < Ef° u /A > is the mean collective energy per nucleon along the direction 



i and q\ is a normalization constant given by [46 



a = t/^f^ [ (6) 

' i=x,y,z 

where N is the number of the particles generated in the simulations. 

It is worth pointing out that the generation of full events is beyond the scope of the 
present simulations. Our purpose here is to simulate in a very simple way the phase space 
distributions of fragments of different masses emitted in the mid-rapidity region. The aim 
is to explore whether the apparent flow angle provided by IMF sideward flow data can be 
considered as an accurate measure of the the effective flow angle. 

The simulations were carried out for semi-central events in the MUL4 event bin at an 
incident energy of E = 150 AMeV where IMFs are copiously produced . We have 

used the following values as input parameters : T = 20 MeV, < E^ oll /A >=< Ey° n /A >= 
3.4 MeV, < E™ u /A >= 8 MeV and 9^ ut = 27.5°. Within this set of parameters, we 
could achieve a quite reasonable reproduction |§ of the dN/p t dp t distributions measured in 
the mid-rapidity region. The azimuthal anisotropy due to the squeeze-out effect [0,0] was 
neglected (i.e. we assumed that < E£ oll /A >=< Ey° u /A >). The value of used in the 

calculations was taken equal to the apparent flow angle measured for Z = 7 fragments. The 
temperature parameter was found to affect only weakly the flow parameter of IMFs (less 
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subject to thermal fluctuations). Its value was constrained to match the experimental flow 
parameter of Z = 2 particles. 

In Fig. 5 we compare the experimental data (open circles), expressed in terms of the 
apparent flow angle versus the fragment charge, to the outcome of the simulations with 
(stars) and without (triangles) taking into account the detector filter [37 1. In the calculations, 
the mass of A > 4 particles was assumed to be equal to twice their charge A = 2 x Z. For 
A < 3 particles, the experimental isotopic ratios f|8| were taken into account. As it can be 
seen, the filtered calculations supply a fairly good quantitative reproduction of the observed 
charge dependence of the apparent flow angle. Note in particular that the saturation seen 
in the data for the heaviest fragments is also present in the simulations. As the particles 
become heavier, the value of the apparent flow angle converges toward the value of the 
effective flow angle (dotted line in Fig. 5) used as an input in the simulations. This 

strengthens the idea that a good estimate of the mean flow angle can be extracted from 
the present measurements of IMF sideward flow. It is also interesting to notice that the 
biases introduced by the FOPI apparatus do not affect significantly the apparent flow angle 
of IMFs. At the incident energy considered here (E = 150 AMeV), above Z = 5 (Fig. 5) 
filtered and unfiltered calculations give almost the same results. 



E. Scaling properties of the flow angle 

It is interesting to investigate now the dependences of the flow angle on the impact 
parameter and the incident energy. To this end, we applied the method described above to 
evaluate the flow angle from the data. Practically, to quantify the constant limiting value 
from the Z— dependence of the apparent flow angle, we calculated the average over the 6^ 
values measured for the heaviest fragments in the saturation region. The range of fragment 
charges included in the average, at each beam energy, is indicated by the shaded area in 
Fig. 3. The results are displayed in Fig. 6 for the Au+Au reaction. Error bars are the 
resulting uncertainties on the average value. 
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As is seen in the left-hand panel of Fig. 6 at an incident energy of E = 400 AMeV, the 
flow angle is found to depend strongly on the collision centrality. We observe a monotonic 
increase of 0fl ow with decreasing impact parameters. This trend is qualitatively in agreement 



with the predictions of several theoretical models p9| , |4*9| . |50|] . The latter models predict 
that Oflow reaches 90° (with an oblate shape), in the limit of the most central collisions 
(b — > 0). This very interesting issue cannot, however, be directly addressed on the basis of the 
present data as the procedure used to evaluate Ofl ow is not reliable for low impact parameter 
collisions. For relatively central events {b geo < 3 fm), the dispersion on the reconstructed 
reaction planes is dramatically large and the effect on the apparent flow angle becomes 



arc 



difficult to be accounted for. Investigations of cm. polar angular distributions |14] , |5T 
expected to be better suited for elucidating this still debated question concerning the flow 
pattern in highly central collisions. 

The excitation function of the sideward flow was already investigated in several experi- 
ments [|],|10],^f40],[52| . However, in all these studies the sideward flow was expressed in terms 



of the flow parameter or the directed transverse momentum (pf r ) which are, in contrast to 
6fl ow , subject to the fluctuations introduced by the thermal motion. Here, we report for 
the first time on the incident energy dependence of the flow angle. The latter observable 
being a measure of the "pure" collective motion, its evolution with the bombarding energy 
is expected to better reflect the collective response of the nuclear system under different 
conditions. The experimental results are presented in the right-hand panel of Fig. 6 for the 
Au+Au reaction over a very broad energy range going from 90 to 800 AMeV. The data are 
shown for semi-central (MUL4 or PMUL4) events where the directed sideward flow is close 
to its maximum. 

The observed trend (Fig. 6) is characterized by a steep rise of 0fl ow to a maximum at 
around E = 250 — 400 AMeV, followed by a moderate decrease as the bombarding energy 
increases further. This is in complete disagreement with the scaling behaviour, i.e. constant 
flow angle, expected from ideal-fluid hydrodynamics P^J49]]. Quantitative calculations in 



the framework of one-fluid viscous hydrodynamics |^9| have shown that Ofl ow scales with the 
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impact parameter almost independently of the incident energy (from 200 to 800 AMeV). 
It was also found in these calculations that the equation-of-state and the viscosity do not 
influence significantly the value of the flow angle. It seems therefore that, within this model, 
the flow angle is mainly governed by the collision geometry. 

Deviations from scaling are visible in the data (Fig. 6) at both low and high incident 
energies. Such deviations are generally thought to signal the presence of phenomena 



known to violate scaling such as the influence of the equation-of-state and possible phase 



transitions pq| . On the low energy end, the departure from scaling is particularly striking. 
One observes a drastic drop of the flow angle below 150 AMeV. This effect was already 
reported in an earlier publication where we have presented the excitation functions of the 
flow parameter ||. It is probably caused by a transition from predominantly attractive to 
repulsive forces resulting in a liquid-to- vapour phase transition. A rough linear extrapolation 
toward lower incident energies leads to an intersection energy where the sideflow vanishes, 
the so-called balance energy, of about 54 AMeV. This value is somewhat lower than our 
previous evaluation in Ref. []3]|§. This is due to the fact that the data in Ref. were 
corrected for the resolution of the reaction plane using the method of Ref. ||. The latter 
method applied at low incident energies E < 120 AMeV, where the fluctuations on the 
reaction plane are relatively large, overestimates the correction factors and leads, therefore, 
to lower flow values |53| . 

Changes in the nuclear viscosity, T], might also contribute to the sudden change of Ofl ow 
at low energies. Below ~ 100 AMeV, r\ is relatively high because of the influence of the 
Pauli principle. With increasing energies, r\ is expected to drop allowing therefore for a fast 



(sudden) buildup of pressure. This was quantitatively investigated in reference |54| where 
realistic viscosities, derived from Uehling-Uhlenbeck equations, were used. 

The origin of the smooth decline of Ofl ow observed at high bombarding energies (Fig. 6) 
is not yet clearly understood. It might be also attributed to a combination of several 
effects : the anisotropy of the nucleon-nucleon interaction, the excitation of resonances and 
possible changes in the stiffness of the equation-of-state |55 |. Recent data from AGS [5£| 



and CERN/SPS H57| , [58f experiments at higher energies indicate relatively small sideward 
flow effects. An estimate of Ofl ow ~ l/20deg. was reported from the Pb(158AGeV)+Pb 



data IB? 



Overall, the rise and fall of the sideward flow, although it is not yet clearly understood, 
reflects certainly important changes in the properties of strongly interacting hot and dense 
hadronic matter. At the low energy end, the onset of flow is indicative of a reduced pressure 
in the system and hence might testify to a liquid-gas phase transition. Similarly, the fall of 
flow on the high energy side might also result from the possible transition towards the quark 



gluon plasma phase which is expected to soften the equation-of-state ||55|| . The present data 
deserve, therefore, to be further investigated in conjunction with the higher energy data from 
AGS and SPS experiments. First attempts along this line have started very recently [ffl. 



V. SUMMARY AND CONCLUSION 

In this work, a systematic study of the directed sideward flow of light charged particles 
and intermediate mass fragments was achieved with the FOPI detector at the SIS/ESR 
facility. Three symmetric reactions, Ru+Ru, Xe+Csl and Au+Au, were investigated at 
beam energies between 90 AMeV and 800 AMeV. The data were analysed according to 
the transverse momentum analysis method and corrections for the finite resolution of the 
reaction plane reconstruction were applied. 

The flow parameter exhibits the expected correlation with the collision centrality with a 
well defined maximum at an impact parameter around b = 4.8 fm and vanishing sideflow at 
both the low and high b ends. Consistent with earlier results [fiOfl , the sideward flow scales 
according to the Ap 3 + A^ 3 scaling rule. 

The flow parameter is found to rise with the charge of the detected fragment up to 
Z = 3 — 4 and then turns into saturation for heavier fragments. This trend supports the 
concept of a collective motion and can be quantitatively accounted for by simple simulations 
based on an anisotropically expanding thermal source. Based on these simulations, we have 
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shown that the value of the flow parameter at saturation provides a good estimate of the 
flow angle 0fl ow m the participant region, offering therefore a new method to evaluate this 
observable in the beam energy range considered here. In contrast to the flow parameter, 
0fl ow has the advantage to be unaffected by the influence of the thermal motion. 

Applying this method, we have explored the dependences of the flow angle on both the 
collision centrality and the bombarding energy. In the impact parameter range (6 > 3 fm) 
where this method is applicable, Ofl ow rises monotonically as the collision becomes more 
central. This trend is qualitatively in agreement with the predictions of several theoretical 
models |19|,^9|,|50| . The excitation function of Ofl ow reveals clear deviations from the flat 
dependence expected from an ideal hydrodynamical scaling. The data exhibit a steep rise 
of ©flow to a maximum at around 250 — 400 AMeV, followed by a moderate decrease as the 
bombarding energy increases further. On the low energy end, the departure from scaling is 
particularly striking. The flow angle drops drastically below 150 AMeV. This effect might 
be caused by a sudden decrease of the internal pressure due to the influence of the attractive 
mean field. A rough extrapolation toward lower incident energies leads to an estimate of 
the balance energy to about 54 AMeV. This value is consistent with the systematics of the 
balance energy ]59fl . The origin of the smooth decrease of 0fl ow observed above 400 AMeV 
is not yet clearly understood. It might be attributed to a combination of several effects : 
the anisotropy of the nucleon-nucleon interaction, the excitation of resonances and possible 
changes in the stiffness of the equation-of-state [[55|| . A quantitative evaluation of these effects 
would require detailed comparisons with the predictions of microscopic transport models. 
Further studies in conjunction with the higher energy data from AGS and SPS experiments 
should also shed some light on the mechanism responsible for the decline of flow at high 
incident energies. 
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TABLES 

TABLE I. Measured cross section (a), geometrical impact parameter (b geo ), reduced geometri- 
cal impact parameter (b g eo/b ma x), accuracy on the determination of the reaction plane (<j(A$r)) 
and the correction factors (< cos(A<I>r) >) for the three reactions under study. The results are 
given for events in the MUL4 bin at an incident energy of E = 400 AMeV. The maximum impact 
parameter, b max , was calculated assuming ro = 1.20 fm. 







Ru+Ru 




Xe+Csl 




Au+Au 


a (barn) 




0.420 






0.523 




0.738 


bg eo (fm) 




2.83 






3.12 




3.61 


bgeo/bmax 




0.26 






0.26 




0.26 


ct(A$ r ) {deg.) 




34.24 






31.09 




22.4 


< cos(A$ R ) > 




0.833 






0.862 




0.901 


TABLE II. Measured cross section (<r), mean geometrical impact parameter (b geo and 
bgeo/b m ax), accuracy on the determination of the reaction plane (ct(A<3?r)) and the correction 
factors (< cos(A<3?r) >) for Au+Au collisions at different beam energies. The results are given for 
the MUL4 (PMUL4 for the 600 and 800 AMeV data) event class. b max was calculated assuming 
ro = 1.2 fm. 


E (AMeV) 


90 


120 


150 


250 


400 


600 


800 


a (barn) 


0.577 


0.720 


0.608 


0.753 


0.738 


0.575 


0.617 


bgeo (fm) 


3.75 


3.98 


3.53 


3.68 


3.61 


3.44 


3.58 


bgeo/bmax 


0.27 


0.28 


0.25 


0.26 


0.26 


0.25 


0.23 


a(A$ R ) (deg.) 


46.2 


39.03 


33.7 


24.7 


22.4 


24.8 


28.8 


< cos(A^r) > 


0.552 


0.758 


0.838 


0.893 


0.901 


0.885 


0.847 
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FIGURES 

FIG. 1. Normalized mean in-plane transverse momentum per nucleon (< >) as a function 
of the normalized rapidity (y^) for particles with charge Z = 1 (full squares), 2 (open triangles), 3 
(dots), 4 (open squares), 5 (full triangles) and 6 (stars). The data are shown for Au+Au (a), Xe+Csl 
(b) and Ru+Ru (c) reactions for the same MUL4 centrality cut and at the same bombarding energy 
of E = 400 AMeV. They are corrected for fluctuations of the reaction plane. For the sake of clarity, 
an offset (const = 0.7 x (Z — 1)/10) is applied to the data points corresponding to a given particle 
of charge Z. Error bars include statistical errors only. 

FIG. 2. Dependence of the normalized flow parameter on the geometrical impact parameter. 
The data are shown at a bombarding energy of 400 AMeV. Error bars are statistical errors 
multiplied by v5?- Geometrical impact parameters are deduced from the measured cross sections 
assuming a sharp-cut-off approximation. 

FIG. 3. Normalized flow parameter versus the charge of the detected fragment. The data are 
shown for the Au+Au reaction at different beam energies, under the MUL4 or PMUL4 (for 600 
and 800 AMeV) centrality cut. Error bars are statistical errors multiplied by \px^ . The solid lines 
are fits to the data using Fermi functions. The shaded area indicates the Z range used to extract 
the average flow angle (see section IV. E) 

FIG. 4. Normalized flow parameter scaled to (Aj/ 3 + A^/ 3 ) versus the charge of the detected 
fragment. The results are shown for the three reactions under study at an incident energy of 
400 AMeV. The MUL4 centrality cut was applied to all 3 reactions. Error bars are statistical 
errors multiplied by \px* . 
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FIG. 5. Apparent flow angle versus the fragment charge for Au(150 AMeV)+Au semi-central 
(MUL4) collisions. Error bars take into account statistical errors only. The data (open circles) are 
compared to the results of anisotropic expanding thermal source model simulations with (stars) 
and without (triangles) the FOPI detector filter. The horizontal dotted line corresponds to the 
value of Ojj^ 1 * = 27.5° used in the calculations. See text for more details. 

FIG. 6. Left-hand panel : Mean flow angle as a function of the geometrical impact parameter. 
Data are shown for the Au+Au reaction at an incident energy of E = 400 AMeV. Right-hand panel 
: Mean flow angle as a function of the bombarding energy. Data are shown for Au+Au collisions 
under the MUL4/PMUL4 centrality cut corresponding to mean geometrical impact parameters in 
the range 3.5 - 4 fm (see Tab. 2). 
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Figure 1 
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Figure 2 
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Figure 4 
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Figure 5 
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Figure 6 
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